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Spherically symmetric oscillatons (also referred to as oscillating soliton stars) i.e. gravitationally 
bound oscillating scalar lumps are considered in theories containing a massive self-interacting real 
scalar field coupled to Einstein's gravity in 1 + D dimensional spacetimes. Oscillations are known 
to decay by emitting scalar radiation with a characteristic time scale which is, however, extremely 
long, it can be comparable even to the lifetime of our universe. In the limit when the central density 
(or amplitude) of the oscillaton tends to zero (small- amplitude limit) a method is introduced to 
compute the transcendentally small amplitude of the outgoing waves. The results are illustrated in 
detail on the simplest case, a single massive free scalar field coupled to gravity. 



I. INTRODUCTION 

Numerical simulations of Seidel and Suenfl| have re- 
vealed that spatially localized, extremely long living, os- 
cillating configurations evolve from quite general initial 
data in the spherically symmetric sector of Einstein's 
gravity coupled to a a free, massive real Klein-Gordon 
field. For example, they observed that initially Gaus- 
sian pulses evolve quickly into configurations which ap- 
pear to be time-periodic. It has been already noted in 
Ref.[l|, that the resulting objects may not be strictly 
time-periodic, rather they may evolve on a secular time 
scale many orders of magnitude longer than the observed 
oscillation period. These interesting objects were first 
baptized " oscillating soliton stars" in Ref. [ij , but some- 
what later the same objects have been referred to as "os- 
cillatons" by the same authors 0. This latter name has 
been by now widely adopted, and we shall also stick to 
its usage throughout this paper. It has been observed 
in the numerical simulations of Ref. [l| that oscillatons 
are stable during the time evolution. Moreover it has 
been argued in Ref. ^ that oscillatons do form in physi- 
cal processes through a dissipationless gravitational cool- 
ing mechanism, making them of great physical impor- 
tance. For example oscillatons would be good candidates 
for dark matter in our Universe. 

On the other hand, stimulated by the seminal work 
of Dashen, Hasslacher and Neveu in the one-dimensional 
(/)^-theory Q, numerical simulations have revealed that 
in an impressive number of scalar field theories spatially 
localized structures -oscillons- form from generic initial 
data which become very closely time periodic, and live for 
very long times [3-[T2|. These objects oscillate nearly pe- 
riodically in time, resembling "true" (i.e. time-periodic) 



breathers. An oscillon possesses a "radiative" tail outside 
of its core region where its energy is leaking continuously 
in form of (scalar) radiation. Therefore a simple approx- 
imate physical picture of a sufficiently small-amplitude 
oscillon is the that of a "true" breather whose frequency 
is increasing on a secular time scale since the amplitude 
of the outgoing radiation is much smaller than that of the 
core. It has been shown in Refs. [l^, [3], that slowly ra- 
diating oscillons can be well described by a special class of 
exactly time-periodic "quasibreathers" (QB). Being time 
periodic, QBs are easier to describe mathematically by 
ordinary Fourier analysis than the long time asymptotics 
of oscillons. A QB possesses a localized core in space 
(just like true breathers) which approximates that of the 
corresponding oscillon very well, but in addition it has a 
standing wave tail whose amplitude is minimized. This 
is a physically motivated condition, which heuristically 
singles out "the" solution approximating a true breather 
as well as possible, for which this amplitude would be 
identically zero. The amplitude of the standing wave tail 
of a QB is closely related to that of the oscillon radiation, 
therefore its computation is of prime interest. Roughly 
speaking "half" of the standing wave tail corresponds to 
incoming radiation from spatial infinity. It is the incom- 
ing radiation that maintains the time periodicity of the 
QBs by compensating the energy loss through the outgo- 
ing waves. In a series of papers [1Jt13 a method has been 
developed to compute the leading part of the exponen- 
tially suppressed tail amplitude of QBs, in a large class 
of scalar theories in various dimensions, in the limit when 
the QB core amplitude is small. Although oscillons con- 
tinuously loose energy through radiation, many of them 
are remarkably stable. The longevity and the ubiquity 
of oscillons make them of potentially great physical in- 
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terest [17h21| . Quite importantly oscillons also appear in 
the course of time evolution when other fields, e.g. vec- 
tor fields are present p2l - [23 |. There is little doubt that 
oscillons and oscillatons are closely related objects. 

The basic physical mechanism for the anti-intuitively 
slow radiation of oscillons is that the lowest frequency 
mode of the scalar field is trapped below the mass thresh- 
old and only the higher frequency modes are coupled to 
the continuum. 

In this paper we generalize the method of Refs. [l^ - 
[Tg} to compute the mass loss of spherically symmet- 
ric oscillations induced by scalar radiation in the limit 
of small oscillaton amplitudes, e, in 1 + D dimensional 
spacetimes. These methods have been succesfully applied 
to D-dimensional scalar field theories coupled to a dila- 
ton field j25i]. Numerous similarities exist between cou- 
pling a theory to a dilaton field and to gravitation: the 
field configurations are of order and the lowest order 
equations determining the profiles are the Schrodinger- 
Newton equations. The stability pattern is also analo- 
gous. Despite these similarities between the dilaton and 
the gravitational theory there are some technical and 
even some conceptual differences. Since there is no time- 
like Killing vector neither for oscillatons nor for the cor- 
responding QBs, already the very definition of mass and 
mass loss is less obvious than in flat spacetime. Another 
conceptual issue is that the spacetime of a time-periodic 
QB is not asymptotically flat, which is related to the fact 
that the "total mass" of a QB is infinite. In the case of 
spherical symmetry considered in this paper a suitable 
local mass function is the Misner- Sharp energy and the 
mass loss can be defined with aid of the Kodama vector. 
The issue of the precise asymptotics of spacetimes can 
be sidestepped in the limit e — > by considering only a 
restricted, approximatively flat spacetime region contain- 
ing the core of the QB (having a size of order 0{l/e)) and 
part of its oscillating tail. We flnd that to leading order 
in the e expansion the oscillaton core is determined by 
the D-dimensional analogues of the Schrodinger-Newton 
equations [26^30'] independently of the self-interaction 
potential. It turns out that exponentially localized os- 
cillatons exist for 2 < D < 6. These flndings show a 
striking similarity to dilaton-scalar theories as found in 
Ref . [2^ . In the case of spherically symmetric oscillatons 
no gravitational radiation is expected due to Birkhoff's 
theorem. The mass loss of spherically symmetric oscilla- 
tons is entirely due to scalar radiation. 

The following simple formula gives the mass loss of a 
small-amplitude oscillaton in D spatial dimensions: 

dM ci / C2\ 

where m denotes the mass of the scalar field, ci is a In- 
dependent constant, while C2 depends on both D and the 
self-interaction scalar potential. The numerical values 
of ci , C2 in the Einstein-Klein-Gordon (EKG) theory for 
spatial dimensions D = 3, 4, 5 are given in Table IVTl We 
also compute and tabulate the most important physical 



properties of oscillatons in the EKG theory (their mass as 
a function of time, their radii). We would like to stress, 
that the method is applicable for oscillatons in scalar 
theories with any self-interaction potential developable 
into power series. 

In the seminal work of Don N. Page [U both the clas- 
sical and quantum decay rate of oscillatons has been con- 
sidered for the case of free massive scalars in the EKG 
theory (for D — 3). We agree with the overall qualita- 
tive picture of the oscillaton's mass loss found in Ref.fsi'], 
however, there are also some differences in the quantita- 
tive results. For example, the amplitude of the outgoing 
wave (related to found by our method differs sig- 

nificantly from that of Ref. f3l| . The main source of this 
discrepancy is due to the fact that this amplitude is given 
by an infinite series in the e expansion, where all terms 
contribute by the same order, whereas in the estimate 
of Ref.fsil only the lowest order term in this series has 
been used. Our methods which are based on the work of 
Segur-Kruskal [13] avoid this difficulty altogether, more- 
over for the class of self-interaction potentials containing 
only even powers of the scalar field, $, the radiation am- 
plitude can be computed analytically using Borel sum- 
mation. 

We now give a lightning review on previous results 
scattered in the literature on oscillatons in 3-1-1 dimen- 
sions. For a given scalar field mass, m, there is a one- 
parameter family of oscillatons, parametrized, for exam- 
ple, by the central amplitude of the field, $c. As $c 
increases from small values, the mass of the oscillaton, 
M, is getting larger, while the radius of the configuration 
decreases. For a critical value of the central amplitude, 
$criti a maximal mass configuration is reached. Oscilla- 
tons with central amplitudes > 'I'crit are unstable 
This behavior is both qualitatively and quantitatively 
very similar to that of boson stars ^2^ i33,, j54.] , and also to 
the behavior of white dwarfs and neutron stars [s^ . For 
reviews of the vast literature on boson stars see for exam- 
ple, Refs. and js^j. In Refs. [H,!!!! a one-parameter 
family of oscillaton-type solutions in an Einstein-scalar 
theory with two massive, real scalar fields has been pre- 
sented, which are essentially transitional states between 
boson stars and oscillatons. 

The interaction of weak gravity axion field oscillatons 
with white dwarfs and neutron stars have been discussed 
in [iol, |4l| , proposing a possible mechanism for gamma 
ray bursts j42l |. Since for very low mass scalar fields os- 
cillatons may be extremely heavy, it has been su gge sted 
that they may be the central object of galaxies [4j|, or 
form the dark matter galactic halos |44-50l. 

Qualitatively good results for various properties of os- 
cillatons has been obtained by Ureha-Lopez [51], truncat- 
ing the Fourier mode decomposition of the field equations 
at as low order as cos(2wt), where w is the fundamental 
frequency. Then the space and time dependence of the 
scalar field separates as ^{t,r) = $i(r) cos(a;i). Oscilla- 
tons with nontrivial self-interaction potentials have also 
been studied in [sij, indicating that similarly to boson 
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stars, the maximal mass can be significantly larger than 
in the Klein-Gordon case. 

The Fourier mode equations have been studied in [s^ 
up to orders cos(10a;i). The obtained value of the max- 
imal mass by this higher order truncation is 0.607/m in 
Planck units. For small-amplitude nearly Minkowskian 
configurations spatial derivatives are also small, and in 
Ref. 52] (and independently in [53]) it has been demon- 
strated that such nearly flat oscillatons can be described 
by a pair of coupled differential equations, the so called 
time independent the Schrodinger-Newton equations [26l - 
|30| . These equations also describe the weak gravity limit 
of boson stars. For quantum mechanical motivations 
leading to the Schrodinger-Newton equations see (54, 55] . 

The time evolution of perturbed oscillatons has been 
investigated in detail by (Sg ■ For each mass smaller than 
the maximal oscillaton mass there are two oscillaton con- 
figurations. The one with the larger radius is a stable S- 
branch oscillaton, and the other is an unstable U-branch 
oscillaton. Moderately perturbed S-branch oscillatons vi- 
brate with a low frequency corresponding to a quasinor- 
mal mode. Perturbed U-branch oscillatons collapse to 
black holes if the perturbation increases their mass, oth- 
erwise they migrate to an S-branch oscillaton. Actually, 
U-branch oscillatons turn out to be the critical solutions 
for type I critical collapse of massive scalar fields (57| . 
Corresponding apparently periodic objects also form in 
the critical collapse of massive vector fields [s^ . 

There are also excited state oscillatons, indexed by the 
nodes of the scalar field. The instability and the decay of 
excited state oscillatons into black holes or S-branch os- 
cillatons is described in [5^ . The evolution of oscillatons 
on a full 3D grid has been also performed in [5£], cal- 
culating the emitted gravitational radiation. Since f{R) 
gravity theories are equivalent to ordinary general rela- 
tivity coupled to a real scalar field, oscillatons naturally 
form in these theories as well [g^] • The geodesies around 
oscillatons has been investigated in [gTI]. 

The plan of the paper is the following. In Section [IT] 
the general formalism concerning a classical real scalar 
field coupled to gravitation in D dimensional, spherically 
symmetric spacetimes is set up. In subsection III CI the 
coupled Einstein-scalar equations are explicited in a spa- 
tially conformally fiat coordinate system. In Section illl] 
the small-amplitude expansion is presented and is carried 
out in detail. In subsection IIII CI it is shown that in lead- 
ing order one obtains the Schrodinger-Newton eqs. in D 
dimensions. In subsection IIII Di thc next to leading order 
results are given. Subsection IIII El contains an analysis 
of the singularities in the complexified radial variable. In 
Section IIVI the proper mass resp. the total mass of the 
QB core is evaluated in subsection IIV Al resp. subsection 
IIVBI In subsection llV DI a conjecture for a criterion of os- 
cillaton stability is formulated. In Section |V] the Fourier 
analysis of the field equations is related to the small- 
amplitude expansion, and the amplitude of the standing 
wave tail of the QB is determined using Borel summa- 
tion techniques. In subsection IV Fl the mass loss rate of 



oscillatons in the EKG theory is computed for _D = 3, 4, 5 
and for various values of the mass of the scalar field. 



II. SCALAR FIELD ON CURVED 
BACKGROUND 

A. Field equations 

We consider a real scalar field $ with a self-interaction 
potential ?7($) in a 13 + 1 dimensional curved spacetime 
with metric gab- We use Planck units with G — c — 
h = I. For a free field with mass m the potential is 
[/($) = m^$^/2. The total Lagrangian density is 



C = Cc 



IGttC 



M ■ 



(2) 



where the Einstein Lagrangian density is Cg — \J~g R, 
and the Lagrangian density belonging to the scalar field 
is 



(3) 



Variation of the action with respect to $ yields the wave 
equation 



(4) 



while variation with respect to yields Einstein equa- 
tions 



Gab — SirTab ■ 
where the stress-energy tensor is 



Tab = $,a$,, 



9ab ( l^.c^^' 



C/($) 



(5) 



(6) 



If Z3 = 1 then, by definition, the Einstein tensor is trace- 
less, and from the trace of the Einstein equations it fol- 
lows that [/(<&) — 0. Hence we assume that D > 1. 

We shall assume that the self-interaction potential, 
U{^), has a minimum C/($) = at <& = 0, and expand 
its derivative as 



(7) 



k=l 



where Uk are constants. In order to get rid of the Sir 
factors in the equations we introduce a rescaled scalar 
field and potential by 



= VStt $ , 



C/(0) = 87r[/($) . 



Then 



fc=i 



(8) 



(9) 
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with 



Vk 



Uk 



(87,)(fc-l)/2 



(10) 



The mass of the field is m = y/ul = ^/vi. If the pair 
4>{x'^) and gab{x'^) solves the field equations with a poten- 
tial C/(0), then = Hl^") and gabix") = gabiix"), 
for any positive constant 7, is a solution with a rescaled 
potential 7^J7(0). It is sufficient to study the problem 
with potentials satisfying = ui — vi — 1, since the 
solutions corresponding to an arbitrary potential can be 
obtained from the solutions with an appropriate poten- 
tial with m = 1 by applying the transformation 

(jyix") (jyimx") , gabix") gabimx") . (11) 

To simplify the expressions, unless explicitly stated, in 
the following we assume m — 1. 



B. Spherically symmetric D -I- 1 dimensional 
spacetime 

We consider a spherically symmetric DA- \ dimensional 
spacetime with coordinates = {t,r,0i, ...,0d-i)- The 
metric can be chosen diagonal with components 



9tt ^ -A , grr^B, 



geiBi = C , 



C Jl sin2 Ok 



(12) 



k=l 



where A, B and C are functions of temporal coordinate 
t and radial coordinate r. The nonvanishing components 
of the Einstein tensor and the form of the wave equation 
are given in Appendix [X] 

A natural radius function, f, can be defined in terms of 
the area of the symmetry spheres in general spherically 
symmetric spacetimes. In the metric (1121) it is simply 



(13) 



C. 

The Kodama vector [6^ [6^ is defined then by 



(14) 



where tab is the volume form in the {t, r) plane. Choosing 
the orientation such that e^t — V AB makes K"' future 
pointing, with nonvanishing components 



K'' = — 



(15) 



It can be checked that, in general, the Kodama vector 
is divergence free, K":^ = 0. Since contracting with the 
Einstein tensor, C^^Ka-b — 0, the current 



Ja = TabK 



(16) 



is also divergence free, = 0, it defines a conserved 
charge. Integrating on a constant t hypersurface with 



a future oriented unit normal vector n", the conserved 
charge is 



E : 



2tt- 
mi 

' D 



r(f)Vo A 



Bn^Jadr 



{Tttf,r - Ttrf^t) dr . 



(17) 



It is possible to show [62, [63], that E agrees with the 
Misner-Sharp energy (or local mass) function rh j64| . 
which can be defined for arbitrary dimensions by 



(£>- l)7r^ 

8-r (f ) 



-D-2 



(1 - g'^'f^^f^b) . (15 



It can be checked by a lengthy calculation, that the 
derivative of the mass function is 



TO , 



2^^f^-i 

' r(f) 



-iabJ 



For the radial derivative follows that 
27rTf^-i 



TO - 



T{§)A 



{Tttr,r - Ttrf^t) 



(19) 



(20) 



which, comparing with (jl7p . gives E = rh. Since for large 
r the function m tends to the total mass, this relation 
will be important when calculating the mass loss rate 
caused by the scalar radiation in Section IV Fl The time 
derivative of the mass function is 



2tt 2 r 



rrif 



T{%)B 



(21) 



This equation is according to the expectation, that, be- 
cause of the spherical symmetry, the mass loss is caused 
only by the outward energy current of the massive scalar 
field. If at large distances the metric becomes asymptot- 
ically Minkowskian, A = B = 1, C = and f = r, then 
using dH) and ©, 



r(f) 



8.r (f ) 



(22) 



C. Spatially conformally flat coordinate system 

The diffeomorphism freedom of the general spherically 
symmetric time-dependent metric form (|12|) can be fixed 
in various ways. The most obvious choice is the use of 
Schwarzschild area coordinates by setting C = r^. How- 
ever, as it was pointed out by Don N. Page in [sij, for 
the oscillaton problem it is more instructive to use the 
spatially conformally flat coordinate system defined by 



C 



r^B. 



(23) 
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even if some expressions are becoming longer by this larger amplitude than the constant r spheres in the con- 
choice. As we will see in Sec. IIIII and in Appendix [B] formally flat coordinate system. In both coordinates, 
inside the oscillaton the spheres described by constant when the functions A and B tend to 1, the spacetime 
Schwarzschild r coordinates are oscillating with much approaches the flat Minkowskian metric. 



J 



In the spatially conformally flat coordinate system the Einstein equations take the form 



(D-l) 



A 



_ 6-D 

B — 



(D-l) 



{D - 2) (r^B) 



r'^AT^B] - 



D -2 



A^B' 



tB / A^r 

A^ KrA^B 



Ai 



D-l .1 f B.t 

^2 



{D-2)rB^ 



A^B 

B,r 

rBi 



A 
B 



('A.r)V|(0,O' 



2AU{4>), (24) 

2BC7(0), (25) 

(26) 
(27) 



The right hand sides are equal to 2Gtt^ 2Grr, Gtr and 2{Geiei/r'^ — Crr), respectively. The wave equation is then 



„2D-2 A 



A 2r'^D-2ABD-i 



B' 



t/'(0)=O. 



(28) 



III. SMALL-AMPLITUDE EXPANSION 

The small-amplitude expansion procedure has been 
applied successfully to describe the core region of one- 
dimensional flat background oscillons in scalar theory 
0, m, [fiHl ■ Later it has been generalized for D -I- I di- 
mensional spherically symmetric systems in p^ . and to 
a scalar-dilaton system in [25]. In this section we gen- 
eralize the method for the case when the scalar field is 
coupled to gravity. 



where e denotes the small-amplitude parameter. We ex- 
pand (j) and the metric functions in powers of e as 

oo 

'/' = E^'''^2fc, (30) 

oo 

A = l + Y^ e'^A^k , (31) 

A.— 1 

oo 

B = l + Y, e^'-Bafc . (32) 



A. Choice of coordinates 

We are looking for spatially localized bounded solu- 
tions of the field equations ([5|) for which cj) is small and 
the metric is close to flat Minkowskian. We use the spa- 
tially conformally flat coordinate system defined by . 
It turns out, that under this approximation, all configura- 
tions that remain bounded as time passes are necessarily 
periodically oscillating in time. We expect that similarly 
to flat background oscillons, the smaller the amplitude 
of an oscillaton is, the larger its spatial extent becomes. 
Numerical simulation of oscillatons clearly support this 
expectation. Therefore, we introduce a new radial coor- 
dinate p by 



Since we intend to use asymptotically Minkowskian co- 
ordinates, where far from the oscillaton t measures the 
proper time and r the radial distances, we look for func- 
tions (/)2fc, A2k and i?2fc that tend to zero when p ^ oo. 
One could initially include odd powers of s into the ex- 
pansions (p0 | -([32 |) . however, it can be shown by the 
method presented below, that the coefhcients of those 
terms necessarily vanish when we are looking for config- 
urations that remain bounded in time. 

The frequency of the oscillaton also depends on its am- 
plitude. Similarly to the flat background case we expect 
that the smaller the amplitude is, the closer the frequency 
becomes to the threshold m = 1. Numerical simulations 
also show this. Hence we introduce a rescaled time coor- 
dinate r by 



p^sr, 



(29) 



T = Lot 



(33) 
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and expand the square of the e dependent factor w as 



k=l 



(34) 



It is possible to allow odd powers of e into the expan- 
sion of w^, but the coefficients of those terms turn out 
to be zero when solving the equations arising from the 
small- amplitude expansion. There is a considerable free- 
dom in choosing different parametrizations of the small- 
amplitude states, changing the actual form of the func- 
tion uj. The physical parameter is not e but the frequency 
of the periodic states that will be given by uj. Similarly 
to the dilaton model in [2^ , we will show, that for spatial 
dimensions 2 < D < 6 the parametrization of the small- 
amplitude states can be fixed by setting w = \/l — e^. 



B. Leading order results 

The field equations we solve are the Einstein equa- 
tions (I24|) - (l27t . together with the wave equation (l28l) . 
using the spatially conformally flat coordinate system 
C = r^B. The results of the corresponding calculations 
in Schwarzschild area coordinates C = are presented 
in Appendix[Bj Since we look for spatially slowly varying 
configurations with an e dependent frequency, we apply 
the e expansion in t and p coordinates. This can be 
achieved by replacing the time and space derivatives as 



d_ 

dt 



d 

OT 



d_ 

dr 



d_ 

dp 



(35) 



and substituting r — p/e. 

From the components of the field equations follows 
that 



i>2 = P2 cos(r + 5) , B2 =b2, 



(36) 



where three new functions, p2, S and &2 are introduced, 
depending only on p. From the e"* part of (1^5)1 it fol- 
lows that (5 is a constant. Then by a shift in the time 
coordinate we set 



5 = 0. 



(37) 



This shows that the scalar field oscillates simultaneously, 
with the same phase at all radii. 

The component of the field equations yield that 



^2 = 02 , 

2 

(f)4 = P4 COS r + ^^^^2. [cos(2r) — 3] 



Bi = b4- 



9 

P2 



A{D-l) 



cos(2r) . 



(38) 
(39) 
(40) 



where 02, P4 and 64 are three new functions of p. If 
D ^2, from the equations also follows that 



a2 



2- D 



(41) 



and that the functions 02 and p2 are determined by the 
coupled differential equations 



02 



dp2 

P2 



D-l da2 
P dp 



D-2 
D- 1 



d2 



D-l dp2 , 

dp'^ p dp 



pI : 



^2) 



(42) 
(43) 



If D = 2 then 02 = 0, and there are no nontrivial local- 
ized regular solutions for &2 and p2, so we assume D > 2 
from now. We note that at all orders sinr terms can be 
absorbed by a small shift in the time coordinate. After 
this, no sin(fcr) terms appear in the expansion, resulting 
in the time reflection symmetry at r = 0. 

Since we have already set rri^ = ui = wi = 1, equations 
(|42l) and (|43|) do not depend on the coefficients Vk of the 
potential U{(j)). To order the functions (j), A and B are 
the same for any potential. This means that the leading 
order small-amplitude behavior of oscillatons is always 
the same as for the Klein-Gordon case. 



C. Schrodinger-Newton equations 



Introducing the functions s and S by 



a;2 - a2 , S = P2 



D-2 
D- 1 



(44) 



equations (|42l) and can be written into the 

form which is called the time-independent Schrodinger- 
Newton (SN) equations in the literature [26. - ,30j : 



d^S D-l dS 



dp^ 



dp 



+ sS = 0, 



d'^s D - 1 ds „n 

— + 5-2 = . 



dp 



p dp 



(45) 
(46) 



Equations P5)) and have the scaling invariance 

{S{p),s{p))^{X'S{Xp),\'s{Xp)). (47) 

If 2 < D < 6 the SN equations have a family of solutions 
with S tending to zero exponentially as p — 00, and s 
tending to a constant sq < as 



So + Sip- 



2~D 



(48) 



The solutions are indexed by the number of nodes of 
S. The nodeless solution corresponds to the lowest en- 
ergy and most stable oscillaton. We use the scaling free- 
dom (j47p to make the nodeless solution unique by setting 
Sq = linip^^oo s — —1. At the same time we change the e 
parametrization by requiring 



CJ2 = -1 for 2 < < 6, 



(49) 



ensuring that the limiting value of 02 vanishes. Then for 
large p 



02 



~sip- 



2-D 



for 2 < D < 6 , 



(50) 
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with only exponentially decaying corrections. Going to 
higher orders, it can be shown that one can always 
make the choice uji — for i > 3, thereby fixing the 
e parametrization, and setting 



= Vf -£2 for 2 < Z) < 6 . 



(51) 



For D = 6 the explicit form of the asymptotically decay- 
ing solutions are known 



s = ±S 



(1 + a2p2) 



for D = 6 , 



(52) 



where a is any constant. In this case, since both s and 
S tend to zero at infinity, we have no method yet to fix 
the value of a in (|52p . Moreover, in order to ensure that 
Lp tends to zero at infinity we have to set 



CJ2 = for L> = 6 . 



(53) 



For D > Q there are no solutions of the SN equations 
representing localized configurations ;66J. 

Motivated by the asymptotic behavior of s, if I? 7^ 2 it 
is useful to introduce the variables 



— , V = s — p a 



2- D dp 



(54) 



In 2 < 13 < 6 dimensions these variables tend exponen- 
tially to the earlier introduced constants 



lim (7 — si 



lim v — sq . 

>-oo 



(55) 



Then the SN equations can be written into the equivalent 
form 



da 

^ 2-D 



52 = 0, 



dv 



dp D-2 



(56) 
(57) 



d^S D-1 dS 



dp"^ 



p dp 



— + U^ + P 



^-°a)S^O, (58) 



which is more appropriate for finding high precision nu- 
merical solutions. Equation ([SG)) will turn out to be use- 
ful when integrating the mass-energy density in Section 
IIV Al in order to determine the proper mass. 



D. Higher order expansion 

From the components of the field equations follows 
the time dependence of A4, 



A4, = af'> + af ^ cos(2t) , 



(59) 



where a^^"^ and a^^'' are functions of p. The functions 
and af"^ are determined by the coupled equations 



dp^ 



D-1 daf'^ _ 2p2Pi{D~2) 



doa 
dp 



P dp 
2 



D- 1 

2pla2 
D-1' 



dp"^ 



D-ldpi 

1— = P4(a2 - ^2) 

P dp 



(^af^ - uji^ p2 



a2P2{D - l){a2 - 1^2) 



D 



(60) 



(61) 



Dpi 



8{D-1) 



5 2 



V3 P2- 



We look for the unique solution for which both a^f^ and 
P4 tend to zero as p — 00. For 2 < D < 6 the function 
P4 goes to zero exponentially, while for large p 



,(0) 



S2P 



S3 : 



(62) 



where si is defined in (|48p. and S2 and S3 are some con- 
stants. If 04"'' and p4 are solutions of (|60l) and (pT|) . then 
for any constant c 



P4=P4 + c[ 2p 



2(a2 - L02] 
dp2 



da2 
dp 



dp 



(63) 
(64) 



are also solutions. This family of solutions is generated 
by the scaling freedom (|Tf|) of the SN equations. If we 
have any solution of (pO)) and (I5T|) then by choosing c 
appropriately we can get another solution for which S3 = 
in 

The equation for 64 is 



db4 
dp 



1 



da 



(0) 



2-D dp 

1 da2 
4(D-2)2 'dp 

P 



da2 
dp 



dp2 
dp 



4(D - 1)02 
2 

-^2(02 - 0J2) 



(65) 



2{D- l){D-2) 
For large p the function 64 tends to zero as 

S2 



6-D 



2 „4-2D 



' 8(13- 2)2 ""i^ ■ 2-D' 
The cos(2r) part of A4 is determined by 



2-D 



dp' 



1 da 



(2) 



P dp 



{D -2){a2-uJ2)pl 
2{D-l) 

D dp2 ( dp2 



(66) 



(67) 



2(D -\) dp \ dp 



D 



-P2 
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We remind the reader, that for 2 < D < 6 the choice 
^2 = — 1, ^^4 = is natural, while for D = 6 necessarily 
UJ2 = 0. For a Klein-Gordon field in Z? = 6 the only 
nonvanishing coefficient is a;4 = — 1. 

Summarizing the results, the scalar field and the metric 
components up to order are 



: e P2 cos T 

jp4 cos T + ^ [cos(2t) - 3] }> + ©(e^) , 



(68) 



4°^ + af ^ cos(2r) +0{e''), (69) 



D-2 

64- 



(70) 



9 

pi 



AiD - 1 



■ cos(2r) 



+ 0(e^). 



Going to higher orders, the expressions get rather com- 
plicated. However, it can be seen that for symmetric po- 
tentials, when V2k — 0, the scalar field (j) contains only 
cos(fcr) components with odd fc, while A and B only con- 
tains even Fourier components. 

Some of the higher order expressions simplifies con- 
siderably when considering symmetric potentials with 
V2k = 0. Because the first radiating mode proportional to 
cos(3t) emerges at order in (/) in symmetric potentials, 
we present its higher order expression for the symmetric 
case 



e p2 cos T + e p4 cos r 

, ( Dpi 
\M{D 



-1 



32 



e>6C0sr (71) 
cos(3r)-|-C'(£8), 



(2)^ 



where is a function of p determined by lengthy differ- 
ential equations arising at higher orders. 

For the Klein-Gordon case in Z? = 3 spatial dimensions 
we plot the numerically obtained functions p2^ 02, P4, 



a 



(0) (2) 



and 64 on Figs. [T]and[2] 



Equations (|68p - ([7ni) determine a one-parameter family 
of oscillating configurations depending on the parameter 
e. This family solves the field equations with a scalar 
field mass m = 1 . By applying the rescaling pT|) to the t 
and r coordinates, we can obtain one-parameter families 
of solutions with any scalar mass m. 

To order, the metric is static. This is the biggest ad- 
vantage of the spatially conformally flat coordinate sys- 
tem C = r^B over the Schwarzschild area coordinates 
C = r^. In the Schwarzschild system the constant r ob- 
servers "feel" an order small oscillation in the metric 
(see Appendix [B|) . The magnitude of the acceleration of 
the constant (r, 0i,02---) observers in the general metric 
(HH) is 



1 dA 



2Ay/B dr ' 



(72) 



which has an order oscillating component when using 
Schwarzschild coordinates, while in spatially conformally 



1.5 




1 1 — 


1 

P2 






D=3 








P4 
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FIG. 1: The exponentially decaying functions p2, pi, and 

(2) 

ajj ' for the small-amplitude expansion of the Klein-Gordon 
oscillaton in the D = 3 case. 




FIG. 2: The functions 02, a^^\ and 64 for the D = 3 Klein- 
Gordon system. These functions tend to zero according to a 
power law for p — >■ 00. 



fiat coordinates the temporal change in the acceleration 
is only of order e^. 

The function W = AB^~^ is equal to 1 to order 
in the conformally flat coordinates. This motivates the 
metric form choice 



ds^ = -Adt^ 



(73) 



— j {dr^ +r^del+r^sm'^9iddl + .. 



which has been employed for the D — 3 case in [3ll |. 



E. Singularities on the complex plane 

As we will see in Section |Vl in order to determine the 
energy loss of oscillatons it is advantageous to extend the 
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functions (f>, A and B to the complex plane. In the small- 
amplitude expansion formalism the extension of the co- 
efficient functions (/)fc , Ak and Bk have symmetrically po- 
sitioned poles along the imaginary axis, induced by the 
poles of the SN equations. We consider the closest pair of 
singularities, located at p — HQd, since these will pro- 
vide the dominant contribution to the energy loss. The 
numerically determined location of the pole for the spa- 
tial dimensions where there is an exponentially localized 
core is 



Qg = 3.97736 , 
Qi ^ 2.30468 , 
Q5 = 1.23595. 



(74) 
(75) 
(76) 



The leading order behavior of the functions near the 
poles can be determined analytically, even if the solution 
of the SN equations is only known numerically on the real 
axis. Let us measure distances from the upper singularity 
by a coordinate R defined as 



p = iQo + R- 



(77) 



Close to the pole we can expand the SN equations, and 
obtain that s and S have the same behavior, 



6 

R^ 



6i{D-l) (D - 1)(L>- 51) 
5QdR 



50QI 



0{R) , 
(78) 

even though they clearly differ on the real axis. We note 
that for D > \ there are logarithmic terms in the ex- 
pansion of s and S, starting with terms proportional to 
R'^lnR. According to (gj) and dM]), the expression ((75)) 
determines the parts of 0, A and B near the pole. 

Substituting into (gOl), jlT]), ([65]) and (HZl), the or- 
der contributions af\ p4, 64 and a^^^ can also be de- 
termined around the pole. We give the results for the 
Klein- Gordon case, when Vk = for > 1: 



,(0) 



9(251? + 208) 32AiD{D - 1) In i? 
' 52{D-2)R^ ^ SbQoiD -2)R3 



a-3 
i?3 



O 



\nR 



P4 



D 



(0) _ 9(43D- 104) 
+ «4 - 26(D - 2)i?4 

HD-\){'dD-S) (l_ 
5Qd{D-2)R3 ^ \R\ 
9(333i:» -I- 832) 324:iD{D - 1) In i? 
260(i5~~2)W ~ SbQoiD - 2)2i?3 
a-3 mi{D-l) 



(79) 



(80) 



(81) 



,(2) 



(D-2)i?3 5Qd{D-2)R^ 

9(6 -L») 
~5(D- 2)i?4 
6i{D - 1){D - 6) 



O 



\nR 



}d{D-2)R^ 



O 



i?2 



(82) 



The constant a_3 can only be determined from the spe- 
cific behavior of the functions on the real axis, namely 
from the requirement of the exponential decay of for 
large real p. 



IV. PROPER AND TOTAL MASS 

A. Proper mass 

In this subsection we present the calculation of the 
proper mass Mp, which is usually obtained by the in- 
tegral of the mass-energy density over a spatial slice of 
the corresponding spacetime. In the next subsection the 
calculation of the total mass M will be performed, by 
investigating the asymptotic behavior of the metric com- 
ponents. The difference Eh — Mp — M defines the gravi- 
tational binding energy, which is expected to be positive. 

The mass-energy density is p — Tabu'^u'', where 
the unit timelike vector has the components 

{1/VA, 0, 0). In terms of the rescaled scalar field (f>, 



1 

8^ 



1 



2A \ dt 



1 



2B \ dr 



(83) 



The total proper mass in the metric (I12p is defined by 
the D dimensional volume integral 



27rt 

r(f) Jo 



drpVBC 



D-l 



(84) 



Applying this for the small-amplitude expansion of oscil- 
latons in spatially conformally fiat coordinates, and using 
that p — er and lu'^ — 1 — e'^, we can write 



Mr, 



2-K^ 



8^r(^i JO 



2(l-l-e2a2) 



£— sm(2r)j +- 



dpp^-\-^{\^eH^- 



2 • , 4 • 

e P2UJ sm T + e P4: sm r 



3 "'F2 

e cos r 



_ dp: 
dp 



(85) 



£^P2 COS T + s'^Pi COS T 



+ £ — — (cos(2t) — 3) 



«2 / 2 \3 

— (e p2 cos t) 



Using (|4T|) and (|43|). for the proper mass we obtain 
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where 



' 8.r(f)7o 
^ 8.r(f)yo 



/ 



(87) 



3D -4 
2{D-2) 



a2P2 



The result turns out to be time independent to this order. 
Although the coefficients Vi of the potential also drop out 
from ((88|) . the dependence on the form of the potential 
still comes in through (j6ip . The leading order behavior, 
(1871) . only depends on the scalar field mass m, which has 
been rescaled to 1 for simplicity. Applying dTTI) to obtain 
solutions with m ^ 1, an rri^ factor appears in (j83p for 
the mass-energy density /i. Since the volume element in 
the integral contains a factor, in all the presented 

proper and total mass formulas an m^~^ factor appears. 

Using , (|55|) and (|56p , the leading order coefficient 

is 



(J - l)7rf 



■ si . 



(89) 



while in the spatially conformally flat coordinate system, 
C — r^B, it can be written as 



4r 



D-2 



'0 



D-2 
''0 . 

4 

1:1-2 



(dr^ + r^dri^j,,) , (91) 



where ro is a constant related to the mass. In general 
spherically symmetric spacetimes it is possible to define 
the natural radius function f by (|13p , and the mass func- 
tion m by p8|) . In both the Schwarzschild and confor- 
mally flat coordinates, for the Schwarzschild- Tangherlini 
metric rh is constant. 



8vrr(f) ° 



(92) 



For the small-amplitude expansion of oscillatons in the 
spatially conformally fiat coordinate system the radius 
function is f = r\/B, where B is expanded according to 
(l32t . Using the rescaled radial coordinate p — er, the 
mass function can be expanded as 



,e~D^^(2) 



0\ 



where 



(93) 



The numerically calculated values for M^^ and Mp'^' for 
the Klein-Gordon field case in various spatial dimensions 
are listed in Table HI 



r(2) 



and 



,(1) 



{D-1)tt-2 



_^dB2 

dp 



(94) 





L» = 3 


D = 4 


D = 5 




1.75266 
-2.11742 
-1.53319 


9.06533 
-43.5347 
-39.0020 


21.7897 
-533.732 
-555.521 



TABLE I: CoefScients of the e expansion of the total mass 
M and proper mass Mp for the m = 1 Klein-Gordon case in 
D = 3, 4, 5 spatial dimensions. 



,(2) 



{D - l)7rf 

8-r(f) 



dB^ 
dp 



dB2 
dp 



D 



The total mass is the limit at 



M = e 



A-D 



4^ dB2 
^ dp 



00, 



(95) 



(96) 



Since i?2 = a2/(2 — D), using the asymptotic form ((50)) 
of a2, we get 



B. Total mass 

Since the scalar field tends to zero exponentially, at 
large distances the metric should approach the static 
Schwarzschild- Tangherlini metric [67| . In Schwarzschild 
area coordinates, with C = r^, this metric has the form 



ds' = -\l- 



„D-2 



,D-2 



dr 



1 



j.D-2 
^-\=2 



dr^ 



r'^dn'^ 



D-l : 



(90) 



■ Sl 



(97) 



agreeing with the leading order coefficient of the proper 
mass, Mp^\ given in ([5^ . Using (|65l) for the derivative 

of 64, 



M(2) = (^-;);; p 



D-l 



da 



(0) 



P 



2{D-2) 



da2 
dp 



dp 



3 da2 



(98) 
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which can be easily calculated numerically, since the ex- 
pression in the limit tends to a constant exponentially. 
The numerical results for the Klein-Gordon field are pre- 
sented in Table ID The proper mass and the total mass 
agree to leading order. However, taking into account the 
next term in the expansion, it turns out that, as it can be 
expected, the gravitational binding energy Ei, = Mp — M 
is positive. 

It is instructive to write the expressions for the total 
mass in natural units, where the oscillaton mass M is 
measured in kilograms while the mass of the scalar field 
m in units of eV/c^: 



M(D=3) = £ (4.66 - 5.63 e^) lO^^fc^ 
(2.94- 14.1 e^) lO^^kg 



M, 



(D=4) 



eV 



eV 



M(c^5) = - (8.63 - 211 e^) lO'^kg 

e V mc 



(99) 
(100) 
(101) 



Since O (e'') terms were dropped, these expressions are 
precise only for small values of e. However, comparing 
to the 3-1-1 dimensional numerical results obtained by 
solving the Fourier mode equations in [52j . it can be in- 
ferred that these total mass expressions give a reasonable 
estimate even when e « 0.5. 



C. Size of oscillatons 

Although oscillatons are exponentially localized, they 
do not have a definite outer surface. A natural defini- 
tion for their size is to take the radius r„ inside which n 
percentage of the mass can be found. It is usual to take, 
for example, n = 95. The mass inside a given radius r 
can be defined either by the integral (|84p replacing the 
upper limit by r, or by taking the local mass function rfi 
in (fT5|) . To leading order in e both definitions give 



M(r) 



(D- l)7r- 



■ cr(er) 



(102) 



where a has been introduced in ([M]) as a function of 
p = er. The rescaled radius p„ can be defined by 



cr(oo) 



n 

100 



(103) 



The numerical values of p„ for various n in D = 3, 4, 5 
dimensions are listed in Table [TTl Restoring the scalar 
field mass m into the expression, the physical radius is 



Pn 

em 



(104) 



In natural units, measuring mc^ in electron volts and r„ 
in meters (Roman m). 



1.97- 10" V 



Pn eV 
e mc^ 



(105) 





D = 3 


D = 4 


D = 5 


Pso 


2.240 


1.778 


1.317 


P90 


3.900 


3.013 


2.284 


P95 


4.471 


3.455 


2.652 


P99 


5.675 


4.410 


3.478 


P99.9 


7.239 


5.692 


4.634 



TABLE II: The radius inside which given percentage of the 
mass is contained for various spatial dimensions. 



Similarly to the total mass expressions in the previous 
subsection, this result is still a reasonable approximation 
for as large e values as 0.5. 



D. Stability 

The stability properties of oscillatons are in many re- 
spects very similar to cold neutron and boson stars. Per- 
fect fluid stars are known to be stable for small pc cen- 
tral densities. As /ic increases, the total mass M also 
increases, until it reaches a maximal value Mmax, where 
according to a theorem in [35|, an unstable radial mode 
sets in. Boson stars have analogous stability properties 
[si]. For oscillatons in the EKG system (for D = 3) 
this behavior has also been observed in Refs. [l| and (56t . 
Oscillatons are closely related to flat background oscil- 
lons, which also behave very similarly to neutron and 
boson stars. Since we use a small-amplitude expansion 
for oscillons and oscillatons it is more instructive to use 
the magnitude of the oscillating central amplitude $c, 
instead of the central density /ic. The central density is 
expected to be a monotonically increasing function of the 
central amplitude. 

Given the very close analogy with oscillons we formu- 
late a general conjecture on the stability of oscillatons. 
We recall that in all known examples the stability pat- 
tern of oscillons is the same, namely if dE/de > os- 
cillons are stable, while when dE/de < oscillons are 
unstable, where E = E{e) is the total energy of the os- 
cillon |14l-[l6l [251 . Therefore we conjecture that the same 
stability pattern holds true for oscillatons, except that 
the energy, E is replaced by the total mass M = M{e) 
of the oscillaton. In other words if the time evolution 
(i.e. energy/mass loss) of an oscillon/oscillaton leads to 
spreading of the core, the oscillon/oscillaton is stable, 
while oscillons/oscillatons are unstable if they have to 
contract with time evolution. Therefore if the conjecture 
is true, the first two terms in the expansion of M enables 
us to determine the stability of oscillatons. 

Taking into account the first two terms in (|96|) , for the 
D = 'd Klein-Gordon field the total mass has a maximum 
at 



-M(i) 



0.525, 



3M(2) 

corresponding to the value of the frequency, w, 



(106) 
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0.851. Although this is just a leading order result, it 
agrees reasonably well with the frequency 0.864 obtained 
by the numerical solution of the Fourier mode equations 
in 15211 . The value of the mass at the maximum is 



-^'^m ax — 



,Af(i) « 0.614. 



(107) 



which is also quite close to the number 0.607 given in 

For an axion with m = 10 ^eV/c^ the maximal mass is 
Afmax — 1.63-10^^A;.g, which is about three times the mass 
of the Earth. The radius of this oscillaton, according to 
the leading order approximation ()105|) . is rgs = 16.8cm, 
while its Schwarzschild radius is 2.42 cm. 

The mass maximum is very important concerning the 
stability of oscillatons. According to [l|,[5^, for the Klein- 
Gordon field in I? = 3 small-amplitude oscillatons with 
£ < Emax are stable, while those with s > Emax are unsta- 
ble, presumably having a single decay mode. For D = 4 
and D — 5 the mass is a monotonically decreasing func- 
tion of £, and all oscillatons are expected to be unstable. 
The situation may be totally different for scalar fields 
with a nontrivial potential U{(p). For example, for V2 = 
and = 4, the coefRcient M^^^ becomes positive, and 
there is no maximum on the energy curve in D = 3 di- 
mensions, consequently even large amplitude oscillatons 
can be expected to be stable. Using this potential in 
D = A dimensions, for small e the mass will be a mono- 
tonically increasing function, so small-amplitude config- 
urations should be stable. For D = 5, in this case, there 
is a minimum in the energy curve, above which stable 
oscillatons can be expected. 



V. RADIATION LAW OF OSCILLATONS 

The methods that we apply in this section for the cal- 
culation of the radiation law of oscillatons have been al- 
ready applied for oscillons formed by scalar fields on fiat 
background. The extension of the Fourier mode equa- 
tions to the complex plane has been first used for the 
one-dimensional 0* theory by Segur and Kruskal [33 |. 
The Borel summation method to calculate the small cor- 
rection near the pole has been introduced by Pomeau, 
Ramani and Grammaticos [g^ . The results has been ex- 
tended to higher dimensional oscillons in fl^ and to a 
scalar-dilaton system in . 



A. Fourier expansion 

Since all terms in the expansion pOI) - ([5^ are exponen- 
tially decaying, the small- amplitude expansion can only 
be applied to the core region of oscillatons. It cannot 
describe the exponentially small radiative tail responsi- 
ble for the energy loss. This is closely related to the 
fact that the expansion is not convergent, it is an asymp- 
totic expansion. Instead of studying a radiating oscil- 



laton configuration with slowly varying frequency, it is 
simpler to consider exactly periodic solutions having a 
relatively large amplitude core and a very small ampli- 
tude standing wave tail. We Fourier expand the scalar 
and the metric components as 



(j) = cos(fcu;t) 



fc=0 



Nf 



A = l + '^Ak cos{kujt) , 

fe=0 
Nf 

B = 1 + ^Bk cos(fcwt) , 



(108) 
(109) 
(110) 



k=0 



where 4>k Ak and Bk only depend on r, and solve the 
Fourier mode equations obtained from Einstein's equa- 
tions and the wave equation. Although, in principle, the 
Fourier truncation order Np should tend to infinity, one 
can expect very good approximation for moderate val- 
ues of Np ■ We assume that the frequency is approaching 
from below the mass threshold m = 1, and in this con- 
text, define the e parameter by e = \/l — w^. 

Regularity at the center require finite values for <j)k , Ak 
and Bk for r = 0, together with 



dr 



= 0, 



dAk 



dr 



0, 



dBk 



dr 



= 0. (Ill) 



Concerning the boundary conditions at ?■ — > oo, it is a 
natural but quite restrictive requirement to assume that 
the metric is asymptotically flat, with the t coordinate 
tending to the proper time for large radii. This implies 
that Ak and Bj. for r — cxo. The Fourier com- 
ponents of the wave equation (|28|) for large r decouple, 
and in this case can be written as 



dr2 



D - 1 d(j)n .,22 IN 

; V in LU - 1) 

r dr 



= 0. 



(112) 



In the relevant frequency range l/2<w<l, ifn> 
2 these equations have oscillatory solutions, behaving 
asymptotically as 



is) 



j.(D-l)/2 
„(=) 



In 



y,(D-l)/2 



sin (rv/n2a;2 - l) (113) 
{r\/~r 



where 7^*^ and 7!^^^ are some constants. According to 



((83l) . this oscillating tail has a mass-energy density, /i, 
proportional to r^^^ . This implies that if ^n'^ or 7!*^^ 
is nonzero for any n > 2, the total proper mass of the 
spacetime is infinite. The requirement of the vanishing of 
all these coefficients together with the central boundary 
conditions are clearly too many conditions to satisfy for 
the given number of second order differential equations. 
In general, regular finite mass exactly periodic solutions 
are not expected to exist. 
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If we require that 7^"-' = and "fn' = for all n 
then all tend to zero exponentially, and we have a 
finite mass asymptotically flat configuration. However, 
in general, this solution will be singular at the center, 
hence we name this solution singular breather (SB). For 
a given frequency, w, the singular breather solution is 
unique by parameter counting. 

Because of their close similarity to oscillatons, it is im- 
portant to study another, presumably unique periodic so- 
lution, the so called quasibreather (QB) solution, which 
is regular at the center, but which has a minimal energy 
density standing wave tail. Our aim is to construct the 
quasibreather solution from the singular breather solu- 
tion. It is important to point out that the quasibreather 
picture is only valid inside some large but finite radius. 
However small the energy density of the oscillating tail is, 
going to very large distances its contribution to the mass 
will not be negligible anymore. Consequently, the as- 
sumption that Ak and Bk tends to zero will not remain 
true for arbitrarily large values of r, and consequently 
Eq. (|112p will also change. For sufRciently large values 
of r, the metric function A increases until it causes to 
change the first radiating mode (either 02 or (f)^) from 
oscillating to exponentially decaying. Increasing r fur- 
ther, all modes will stop oscillating one by one. This 
way we obtain the exactly time-periodic but infinite mass 
"breathers" which are described in details in Section IX 
of the paper of Don N. Page [31]. The quasibreather 
can be considered to be the part of such an infinite mass 
" breather" containing the core and a large portion of the 
tail where the first radiating mode oscillates, requiring 
that the mass inside this region is dominated by that of 
the oscillon core. Since the core amplitude is of the or- 
der e^, while the tail is exponentially suppressed in e, 
the quasibreather picture is valid in a sufficiently large 
volume. 

Since the amplitude of the oscillating tail of the quasi- 
breather is very small, apart from a small region around 
the center the core of the QB is very close to the cor- 
responding singular breather solution. In particular, the 
SB and the QB solutions have the same e expansions. 
For the SB solution the small-amplitude expansion will 
not be valid in a region near the center r = 0, while 
for the QB solution it will fail for large radii where the 
oscillating tail becomes dominant. The size of the re- 
gion, (0,rciiff), around ?- = where the difference be- 
tween the SB and the regular core becomes relevant is 
TdiS = 0{e~^/^) (with S being a constant), whereas the 
size of the SB or QB core is proportional to 1/e. Outside 
of this region, i.e. for r > r^is the difference between the 
SB and the QB will be very small since the singular mode 
turns out to be proportional to the tail amplitude, which 
is exponentially small in terms of the small parameter e, 
while the core amplitude is of order s'^. 

For potentials U{(t>) which are symmetric around their 
minima, i.e. V2k = for integer k, the Fourier expansion 
of the scalar contains only odd, while that of the metric 



components only even terms, 



02k 



, A2k+i = , 



B 



2k+l 



= 0. 



(114) 



For symmetric potentials the first radiating mode is ^3. 
In this section we will concentrate mainly on the Klein- 
Gordon scalar field with = for fc > 1. It is straight- 
forward to generalize the results for symmetric potentials. 

For small-amplitude quasibreather or singular breather 
configurations we can establish the connection be- 
tween the Fourier expansion (|108[) - (|110p and the small- 
ampfitude expansion by comparing to (l55|) -d7T |) . 

For symmetric potentials we obtain: 



( Dpi 

y64:{D~l) ' 32 



(115) 



£>3 = £ 



3 (2) 
'2 + + P2Q4 



+ 0(e'), 



A2-eM'^+0(£^), 



B2 



.4 pI 



4iD-r 



+ 0(£^). 



(116) 
(117) 
(118) 
(119) 

(120) 



B. Expansion near the pole 



As £ — > the amplitude of all Fourier coefficients tend 
to zero. However, extending them to the complex plane, 
for small e they all have pole singularities on the imagi- 
nary axis at r = zLiQu/s, corresponding to the poles of 
the Schrodinger-Newton equations at p = zLiQu, as it 
was discussed in Sec. IIII El As e tends to zero, the poles 
move further and further away from the real axis, but 
close to them the Fourier components 4>k , Ak and Bk are 
not getting small, in fact they have e independent parts. 
We introduce a shifted radial coordinate, y, for an "inner 
region" around the upper pole by 



iQi 



(121) 



The coordinate y is related to the one (R) defined in 
Eq. ([77)1 hy R — ey. Substituting the small-amplitude 
expansion resuhs ([THll-dMl) into (|TT5]) - P^ . and taking 
the limit e -> 0, it follows that in the Klein-Gordon case. 
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near the upper pole 



6 



9991? 



D- 1 
£1-2 



^0 

Bo 

B2 



y2 52{D ~ 2)y4 

27{7D- 12) 
A{){D-2)y^ ^ 
6 9(25i:» + 208) 
f~ 52(Z)-2)y4 

9(6 - 
"5(Z?-2)2;4 +■■•' 

6 9(333L' + 832) 

- 2)?;2 + 260(£i - 2)2^4 
9 



D - 1 
iJ~~2 



(122) 

(123) 
(124) 
(125) 
(126) 
(127) 



We note that since (|122p - (|127l) are expansions in l/y^, 
they are vahd for large y values. In contrast, (|78p - ([5^ 
were calculated assuming small R. Both of these condi- 
tions can hold simultaneously, since R — ey. 

Expressions ()122|) - (I127|) can also be obtained by look- 
ing for the solution of the Fourier mode equations in the 
e — limit near the pole as a power series expansion in 

l/2/^ 



oo ^ 

j=k+i 

oo 

E 

oo ^ 



2k- ^ "2fc-2j 
j=k+l ^ 



(128) 
(129) 
(130) 



^2fe 



]=k+l 



and are constants. The mode 



(i) 

where ip2k+i^ 

equations that we have to solve can be obtained from the 
Einstein equations ([M)) - (l?7l) and from the wave equation 
(Pg]) by substituting pIIg)) -(ITTU l) . Equations (I111)-(I2H1) arc 
not independent. The wave equation follows from the 
Einstein equations by the contracted Bianchi identities, 
and the (t, r) component (j26p is a constraint. The trunca- 
tion of the Fourier expansion at a finite Np order makes 
the mode equations mutually contradictory. However, if 
we choose any three field equations from the 
arising mode equations will clearly have solutions. We 
have checked that our results for the energy loss rate of 
oscillatons are the same for different choices of the three 
field equations. We have also tested that the violation 
of the mode equations obtained from the other two field 
equation tends to zero quickly as Np increases. 

Substituting into the field equations (p4t - ((28)) 

and taking the e — ?> limit close to the pole, some terms 
with lower powers of r can be neglected. Then, inserting 
the l/?/2 expansion (|128p - (|130p into the resulting mode 
equations, because of the omission of odd powers of 1/y, 
the only ambiguity arises at the choice of the signature 
ofip^^\ The calculation of ([T^ - p??)) using the Fourier 



mode equations is technically more simple than using the 
small-amplitude expansion method, and can be done by 
algebraic manipulation programs to quite high orders in 
1/2/- 

Apart from an overall factor, the leading order behav- 
ior of the coefficients V'i"'' j and Z?^"' for large n can 
be obtained by studying the structure of the mode equa- 

(n) 

tions. It turns out that for large n, '03 dominates among 
the coefficients. For the third Fourier mode of the Klein- 
Gordon field. 



kni-iy 



pn-l) 
8" 



3(91?- 10) 
2{D - 2)n 



3(91?- 10) (7L'- 



1{D 



(131) 



where kp, is a factor depending on D and Np. All other 
coefficients grow slower with n asymptotically. Although 
the 1/n and l/n? correction terms may depend on the 
choice of the scalar potential, the leading order behavior 
is the same as in (jl3ip for any symmetric potential. The 
value of the constant kp) will turn out to be crucial for 
the determination of the energy loss rate of oscillatons. 
Calculating the coefficients up to order n = 100 and tak- 
ing into account Fourier modes up to order Np = 6, in 
the Klein-Gordon case we obtain 



fca = -0.301, 
ki = -0.134, 
fc5 = -0.0839. 



(132) 
(133) 
(134) 



C. The singular breather solution near the pole 



Expansion (|128p - (|130p gives an asymptotic series rep- 
resentation of the Fourier components (pk, Ak and Bk. 
The results (|122l) - (|127p can be considered as boundary 
conditions for the Fourier mode equations for 



\y\^ oo , -tt/2 < arg ?/ < 0, 



(135) 



ensuring a unique solution for the "inner problem" . This 
corresponds to the requirement that cf) decays to zero 
without any oscillating tail for r ^ oo along the positive 
half of the real axis, i.e. we consider a singular breather 
solution. 

The Fourier components of the wave equation (j28p can 
be written as 



(P(j)n D -I d<f)n .,22 

dr'^ r dr 



(136) 



where Fn contain nonlinear polynomial terms in 0^ , A^ , 
Bk and their derivatives for k < Np. Using the y coor- 
dinate near the pole and taking the e — > limit. 



+ (n^ - 1)0„ 



(137) 
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where F„ denotes the e — ?> hmit of F„. On the imagi- 
nary axis the 1 jy^ expansion gives real valued functions 
to all orders. As singular breather solutions of the mode 
equations, with boundary conditions (|122p - (|127p in the 
region given by (|135p . the functions (\)^ can have small 
imaginary parts on the imaginary axis, satisfying the left 
hand side of (jl37l) to a good approximation. For sym- 
metric potentials the first radiating component is (^3 . The 
singular breather solution can have an exponentially de- 
caying small imaginary part on the imaginary axis, 

Im 0f ^ = 1/3 exp(-iN/8?/) for Re y = , (138) 

where 1^3 is some constant. On the other hand, since 
the quasibreather solution of the mode equations is reg- 
ular and symmetric, has zero imaginary part on the 
imaginary axis. 

For symmetric potentials the value of 1^3 can be ob- 
tained iDy Borel summing the series (|128p for ^3 [69| . 
The first step is to define a Borel transformed series by 

00 , (n) 

^W-E&i-'"- (139) 

n=2 ^ >' 

Then the Laplace transform of Viz) will give us the Borel 
summed series of ^ which we denote by ^ (y) , 

'^f''(y) = ^°°dte-vQ) . (140) 

We are only interested in the imaginary part of 0f^(2/) 
on the negative imaginary axis, y — —iyi, where yi > 
real. Then the argument of V is z — t/y = it/yi, which 
is pure imaginary with positive imaginary part. Since all 
terms in p39p contain even powers of z, no individual 
term gives a contribution to Im ^ on the imaginary 
axis. The value of Im^f^ is determined there by the 
leading order large n behavior of the series (|139l) . Using 
P3ip and including a term proportional to z^, 

(141) 

where the sign denotes equality up to terms that do 
not give contribution to the imaginary part of (f>f^ on 
the imaginary axis. Transforming the argument of the 
logarithm into product form, only one of the factors gives 
a contribution. 

For purely imaginary y. 

In this case, for t > yiV& we have to integrate along the 
branch cut of the logarithm function. In order to see 



how to go around the singularity at t = yi\/8 we note 
that according to ([^5]) . the 1/y^ expansion P^ - p^ 
has been applied for y — yr — Wi, where y^ and yi are 
positive and real. This corresponds to the requirement of 
exponential decay for r > along the real r axis. Then 

iz = * ^ {-yi + iyr) , (144) 

Vr + Vi 

which shows that the argument of the logarithm in (|142p 
has to go around the singularity in the upper half of the 
complex plane. This means that we approach the branch 
cut of the logarithm at the negative part of the real axis 
from above, where its imaginary part is tt. Then for 
purely imaginary y we can evaluate the imaginary part 
of the integral (|140p by integrating on the branch cut, 

Im^f^(y) = - r dte-*^ 

= -^exp(-^V8y). (145) 

The logarithmic singularity of V (tjy) does not con- 
tribute to the integral. Comparing with (|138l) . 

V3 = -\kD^- (146) 

For asymmetric potentials the leading order radiating 
component will be in ^2, and 

Im = vi exp(-i%/3y) for Re y = . (147) 

Because the dominant behavior of c^g j it is not possible to 
determine the constant vi by the Borel summation. Its 
value can be calculated by numerical integration of the 
Fourier mode equations, following the method presented 
in [32i] and tl5,]. 

It is reassuring that even though we work with a trun- 
cated set of mode equations Birkhoff 's theorem still holds 
in the following sense. Neither a^,"-* nor /J^"-* has an ap- 
propriately singular behavior so that they generate an 
imaginary correction on the imaginary axis for and 
Sfe. The Borel summation procedure does not produce 
gravitational radiation. 

D. Construction of the quasibreather 

As we have already discussed, both the singular 
breather (SB), and the quasibreather (QB) solutions 
are well approximated in a large domain by the small- 
amplitude expansion. Since the tail is exponentially sup- 
pressed in e, apart from a small central region around 
r = 0, where the SB solution gets too large, the QB 
and SB solutions are extremely close to each other. We 
denote the difference in the first radiating Fourier com- 
ponent </)3 of the two solutions by 

c = - 4>r ■ (148) 
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The small function 0™ solves the linearization of the wave 
equation around the singular breather solution. To lead- 
ing order in e this reduces to the flat background wave 
equation 



D - 1 



= 0. 



dr^ r dr 

The general solution of (|149p can be written as 

in, 



(149) 



^_D/2-l 



(150) 

where J and Y are Bessel functions of the first and sec- 
ond kinds, and /3d are constants. The asymptotic 
behavior of the Bessel functions is 



2 ( VK 

cos \x — 

TTX V 2 



sin 



(151) 
(152) 



for arga: < tt and |a::| — >■ oo. The constants au and /3d 
describe the amplitude of the standing wave tails in 01", 
since for large distances from the center, 



1 



r(D-i)/2 



ajj sm 



(153) 



-;3dcos VSr -{D- 1)^ 



Since the SB solution is exponentially decaying, this will 
also be the tail of the QB configuration. 

The second term in (I150p gives a purely real contribu- 
tion to ij)^ on the imaginary axis. However, converting 
(|152p into exponential form. 



1 



exp 



^(2. + 3) 



-IX + — (2t/ + 3) 



(154) 



Substituting the value of 1^3 from (|146l) , obtained by the 
Borel summation, 



fcnTT 



Qd 



exp 



D 



(157) 



For any value of Pd the second term in (|150p gives 
a regular symmetric contribution to 0™, which does not 
change the behavior of the imaginary part on the real 
axis. However, as it is apparent from (I153p . any nonzero 
necessarily increases the tail amplitude, and conse- 
quently the energy density in the tail as well. Hence, in 
order to obtain the minimal tail quasibreather, we set 



/3d = 0. 



(158) 



The standing wave tail of the quasibreather in the 
asymptotic region is given by the first term of (|150p . 

0«^-^V^^^yD/2-i(%/8r)cos(3r) (159) 



^(D-l)/2 



sm 



\/8r- (£>- 1)^ cos(3t). 



Subtracting the regular solution involving the Bessel 
function with a phase shift in time, we cancel the 
incoming radiating component, and obtain the radiative 
tail of the oscillaton, 



«D 



.D/2-1 



Y, 



D/2- 



_i(%/8r) cos(3r) 



- ^D/2-i(V8r) sin(3T) 



a, 



VSr - {D - 1)^ - 3t 



(160) 



Equations ()159p and (|160|) are valid for symmetric po- 
tentials. In both cases, the amplitude of the tail of (p 
at large r is given by a^. According to ([8]), the physi- 
cal amplitude is ajj/^/Sn. Since the transformation (jlip 
changes the coordinates, scales as m*^^"^-'/^ with the 
scalar field mass m. 



we see that the first term in ()150p yields an exponen- 
tially behaving imaginary part along the imaginary axis. 
Close to the upper pole at iQo/s, using the coordinate 
y defined in (|121l) . to leading order in e we obtain that 



Im0- = 



an 



e 



exp 



V8Q1 



-tVSy], (155) 



for Re y = 0. Since our aim is to obtain a non-singular 
QB solution which is symmetric for r — > — r for real r, 
(I155P must cancel the exponential behavior of lm(f>f^ 
given by (|138p . This fixes the amplitude ajj, 

}d\^ _ ( \/8Qd\ 



-2i.3 



exp 



(156) 



E. Tail amplitude 

The scalar field tail calculated in the previous subsec- 
tion is so small that it is not surprising that it has not 
been detected by numerically solving the Fourier mode 
equations in [l| and [l^- In order to relate the magni- 
tude of the oscillating tail to the central amplitude, we 
represent 4> in the core region by <^ = £^P2 cos r, and in 
the tail by (I159P . The tail starts to dominate at a radius 
r — Tt where 



0,r) = e^p2{er) = e^S{er) 



D - 1 



D 



(161) 



equals to aj-,A^ Since s ~ —l + sifp' ^ for large p, 

the asymptotic behavior of S in the relevant dimensions 
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siisi-2) 
8p 



Sd=4{p) = St 



,3/2 



P 

Sd=5{p) — St — 5- 
P 



(162) 
(163) 

(164) 



where the constants si and St are given in Table Hill The 





D = 3 


D = 4 


D = 5 


Sl 

St 


3.505 
3.495 


7.695 
88.24 


10.40 
23.39 



TABLE III: The numerical values of the constants si and St 
in 3, 4 and 5 spatial dimensions. 

values of rt and the amplitude of the tail at that radius, 



"D't 



„(l-£')/2 



/vStt for several e are given in Table 



e 




D = 


3 






D = 


4 






D = 


5 




rt 


$t 


rt 


^-t 


rt 


$t 


0.1 


1160 


8.96 


• 10" 


52 


648 


2.76 


10" 


32 


346 


4.42 


10" 


20 


0.2 


302 


4.63 


■ 10" 


27 


168 


1.06 


10" 


17 


92.6 


6.01 


10" 


12 


0.3 


140 


9.28 


• 10" 


19 


78.2 


9.45 


10" 


13 


45.0 


3.83 


• 10" 


-9 


0.4 


81.6 


1.40 


■ 10" 


14 


46.4 


3.07 


10" 


10 


27.9 


1.03 


• 10" 


-7 


0.5 


54.3 


4.68 


■ 10" 


12 


31.5 


1.03 


■ 10 


-8 


19.8 


7.56 


• 10" 


-7 


0.6 


39.2 


2.30 


• 10" 


10 


23.2 


1.09 


• 10 


-7 


15.1 


2.87 


• 10" 


-6 


0.7 


29.9 


3.76 


• 10" 


-9 


18.0 


5.91 


■ 10 


-7 


12.2 


7.42 


• 10" 


-6 


0.8 


23.8 


3.08 


• 10- 


-S 


14.6 


2.12 


• 10 


-6 


10.3 


1.51 


• 10" 


-5 



TABLE IV: The radius rt where the oscillating tail starts to 
dominate, and its amplitude $t there. 

The tail amplitude should be compared to the central 
amplitude 



1 



Sc 



D-2 



(165) 



where the constants Sc and $ic are given in Table |Vl 
The radius rt where the tail starts to dominate is much 





D = 3 


L> = 4 


D = 5 


Sc 


1.021 


3.542 


14.02 


$lc 


0.288 


0.865 


3.229 


Ph 


2.218 


1.357 


0.763 



TABLE V: The numerical values of the constants Sc, $ic 
and ph, which determine the central amplitude $c and the 
characteristic size rn- 

larger than the characteristic radius of the core, which 



can be defined as the radius r/j where <& = $c/2. Clearly, 
rh = Ph/£, where ph is the value of p for which S = Sc/2. 
The value of ph for various spatial dimensions D is also 
given in Table [V] Clearly, there is only some chance 
to numerically observe the tail for as large e values as 
0.5, which for Z? = 3 is close to the maximum value 
Emax ~ 0.525. It can also be observed from Table IIVI 
that for larger spatial dimensions the radiation is signif- 
icantly stronger. It would be reasonable to first make 
the numerical analysis for = 5, since then the tail has 
the largest amplitude. Even if the Klein- Gordon oscilla- 
tons are unstable in I? = 5, as we have seen in Section 
IIVBI there are scalar potentials, for which large ampli- 
tude oscillatons are stable. Since the exponent in (|157p is 
potential independent, in general, we expect to get a tail 
amplitude of similar magnitude as for the Klein-Gordon 
field. 



F. Mass loss rate 

Since at large distances from the center the mass func- 
tion m agrees with the total mass M, the mass change 
rate of the oscillaton can be calculated from the energy 
current carried by the wave (|160p using (P^ . Averaging 
for an oscillation period. 



dM 



■ exp 



C2 



(166) 



where the D dependent constants are 



Cl 



4r(f) ' 



C2^2V8Qd. (167) 



The numerical values of ci and C2 for various spatial di- 
mensions are listed in Table I VI I The values of C2 are the 





D = 3 


D = 4 


D = 5 


Cl 
C2 


30.0 
22.4993 


7.23 
13.0372 


0.720 
6.99159 



TABLE VI: The constants ci and C2 in the mass loss rate 
expression p66[) for _D = 3, 4, 5 spatial dimensions. 

same for any symmetric potential, but the numbers given 
for Cl are valid only for the Klein-Gordon field. 

The higher e is, the more chance we have to observe the 
presumably tiny energy loss. As we have seen in subsec- 
tion HVBl for D = 3 spatial dimensions oscillons are sta- 
ble for e < emax ~ 0.525. The total mass is maximal at 
£max- Restoring the scalar field mass m into the expres- 
sions, the maximal mass value is Mmax = 0.614/to. Sub- 
stituting into (jl66l) . for the maximal mass Klein-Gordon 
oscillaton we get 



1 dAf 
M "dT 



= -4.3- 10"^ V. 



(168) 
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This expression is valid in Planck units. Expressing M 
in kilograms, and mc^ in electron volts the maximal os- 
cillaton mass is 



M„,ax = 0.614 TOpSp/m = 1.63 • W^kg 



eV 



(169) 



where the Planck mass is mp = 2.18 • lQ~^kg and the 
Planck energy is Ep — 1.22 • IQ^^eV . Expressing t in 
seconds, since tp — 5.39 • lO^'^'^s, for the maximal mass 
configuration we get 



1 AM 
M "dT 



0.066 mc^ 



eV 



(170) 



Note that we only determined the leading order result for 
the radiation amplitude. We saw hints that the small- 
amplitude results give sensible answers for moderate e 
values, however we do not have a good control over non- 
leading terms in the radiation amplitude. In general for 
such high e values we expect to have the same ex pon ential 
factor, but with a different prefactor ci [TBI. [l6ll25j. The 
order of magnitude should be nevertheless correct. 

According to ((96|) . for D ~ 3 spatial dimensions, to 
leading order the total mass is proportional to the am- 
plitude. M = eM^^ym, where from Table [H M^^l = 
1.75266. Substituting into (jl66p . for the Klein-Gordon 
case this yields 



dM 
"dT 



C3 



exp 



\ uiMJ 



where 



C3 = 92.2 



C4 = 39.4337 . 



(171) 



(172) 



Expression (11711) has the same form as the classical mass 
loss formula (122) of 3l|, although the constant corre- 
sponding to C3 is much larger there, it is 3797437.776. 
This means that the amplitude of the radiating tail of 
the scalar field $ is overestimated by a factor 202.9 in 
[3l| . In order to understand the reason for this large 
difference, and why it is necessary to use our more com- 
plicated approach to obtain a correct mass loss rate, we 
first describe the method of [3l| in our formalism. 

For a D = 3 Klein- Gordon system let us consider the 
third Fourier component (|136p of the wave equation ([25)1 
in the a; 1 hmit. Taking the results (|115l) - (|120l) of 
the small-amplitude expansion, and substituting into the 
nonlinear terms on right hand side of (|136p . we obtain 
an inhomogeneous linear differential equation for 03, 



-FT + - + 8(/>3 = P{r) ■ 
ar^ r dr 



(173) 



Here the function P{r) is given by the small-amplitude 
expansion, in a power series form in £, 



Pir) ^yP2k{r)s 



2k 



k=3 



(174) 



For the D 
term is 



3 Klein- Gordon system the leading order 
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(2) 



(er). 



(175) 



Equation p73)) with P(r) = e^Peir\_but setting af ^ = 
0, corresponds to equation (52) of [3l|. Since there the 
Fourier modes are defined in terms of exponentials in- 
stead of cosine functions, the coefficient of the term 



is 3/4 in [3l| instead of 3/16. The term containing 



(2) 

a\ is missing there because of the assumption that the 
gtt = ~A metric component is time independent. How- 
ever, at order either gtt becomes oscillatory, or the 
spatial metric ceases to be conformally flat. 

The oscillating tail responsible for the radiation loss in 
the (f>3 mode can be estimated by integrating (|173p using 
the Green function method, 

Mr) = ^"'^^''^ r fsm{V8f)P{f)df (176) 
V8r Jo 



3in(V8r 



r cos 



The oscillaton core is exponentially localized, hence in 
the tail region it is a very good approximation to write 



03 (r) = a 



cos(-\/8r) 



(177) 



where the constant determining the amplitude is 
1 

a^—j rsin{VSr)P{r)dr . (178) 
v8 Jo 

Since the functions describing the oscillaton by the e ex- 
pansion are symmetric around r = 0, and since the sine 
function can be written as the difference of two exponen- 
tials, 



1 



2iVS 



rexp{V8ir)P{r)dr . (179) 



The functions in the small-amplitude expansion depend 
directly on the rescaled radial coordinate p = er, so it is 
natural to write the integral into the form 



1 



a 



J- 



pexp 



•ip 



p 



[^)dp. (180) 



This can be replaced by a contour integral around the 
upper plane, and can be approximated by taking into 
account the pole which is closest to the real axis. The 
position of the closest pole is the same as that of the 
SN equations p5t and (|46|) . it is at p = iQ^ on the 
imaginary axis. Let us first calculate the contribution 
from the leading Pq term, given by (|175l) . Then, since 
for three spatial dimensions P2 — \/2S, using ([78| and 
I, the behavior near the pole is 



P 



e'Pe 



5i?6 



(181) 
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where R = iQs — p. The omission of the term containing 
from (|175p results in a value which is 5 /3 times that 



-,(2) 



of (|18ip . The residue can be calculated by integrating 
by parts five times, 



233''\/27rQ3 
Pi 



■ exp 



(182) 



This amplitude is much larger than the amplitude 
calculated by the Borel summation method in (|157p . 



a 
as 



2%V2 



121.8. 



(183) 



If we omit the term containing 04 from (|175[) , we have 
a 5/3 factor, and get a/a^ = 202.9, which is the ratio 
of the tail amplitude of 3l| to our value, as we have 
mentioned after Eq. (|172p . 

The fundamental problem with the above calculated 
tail amplitude a is that it is just the first term of an infi- 
nite series, of which all terms give contributions which are 
the same order in s. This can be illustrated by calculating 
the contribution of the next term in P{r). Although the 
small-amplitude expansion yields a rather complicated 
expression for ^8(^)1 it contains terms proportional to 

P2I2, 75404^^ and P20'2o!i \ which have eighth order poles. 
When calculating the integral (|180p one has to integrate 
by parts seven times, so the result will have the same e 
order as the earlier result calculated from the Pg (r) term. 
Even if we could calculate higher order contributions, we 
have no reason to expect that the series converges, and 
even if it would be convergent it may not give a correct 
result for the mass loss. This has already been demon- 
strated for the simpler system of a real scalar field with a 
nontrivial interaction potential on flat Minkowski back- 
ground. It was first pointed out in [70| that there are 
too many boundary conditions to satisfy when solving 
the Fourier mode equations in order to find periodic lo- 
calized breather solutions. In the energy loss rate 
of the long living oscillon configurations was estimated 
by a method analogous to that of j3l|. However, cal- 
culating higher order contributions, it turned out that 
the method gives an incorrect nonzero result even for the 
periodic sine-Gordon breather. Moreover, the expansion 
is not convergent for the 0"* scalar theory. The proper 
approach to calculate the lifetime of oscillons has been 
worked out by [32] and , using complex extension and 
Borel summation. As a result of the above arguments, 
the correct value of C3 in the mass loss rate expression 
([TTTI) for the £> = 3 Klein-Gordon field is cg = 92.2. 

Although the expression (|17ip is correct for small val- 
ues of M, since it is based on the assumption that M 
depends linearly on e, one should not apply it to mass 
values close to Mmax- For example, substituting the value 
of the maximal mass Afmax — 0.614/m into (|17ip . we ob- 
tain 



' 1 AM 
M "dT 



= -5.0- 10^26^, 



(184) 



which is 9 magnitudes smaller than the maximal mass 
loss rate obtained in (|168p . The reason for this huge 
difference is that according to the linear expression M = 
£M'^^\ to the mass value Afmax — 0.614/m belongs an e 
value of 0.350. At that e we obviously get a significantly 
lower radiation than at Cmax ~ 0.525, because of the 
exponential dependence. Since expression ()168p does not 
involve this approximation, we expect it to give a more 
reliable result. 



G. Time dependence 



Instead of using (|171l) to determine the time depen- 
dence of the oscillaton mass, in order to obtain results 
that are valid for larger mass values, we work out a 
method involving a higher order approximation for the 
£ dependence of the mass. Since the first two terms of 
(|96p determine the mass maximum to a good precision, 
we expect it to be a reasonable approximation for close 
to maximal e values. Including the scalar field mass m, 
we use 

Taking the time derivative and comparing with (|166p . 

^ = --(/3i+/32£^)exp(^), (186) 

where 

/3i ^ i— ^M(i) , /32 = ^m(2) . (187) 



Cl 



Cl 



This can be integrated in terms of the exponential inte- 
gral function. 



t-tn = -- 



120to 

- 132 (4 + 



[20/3i {cl + C2£ + 2e^ 



cle- 



2cl£^ + 6c2e^ 

2\ T^. /C2 



(20/3i + p2cl) Ei ^ 



120m ^ -' \ e 

Taking the expansion of the result, for small e 



24e4)] exp (^) 
(188) 
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m 



C2 



+ 



20^1 + 1324 



6e^ 

C2 



exp 



(189) 

(f)- 



M=M„, 



Although the correction from the subleading term M'^^) 
only appears in the third term in the bracket, its influence 
for e « 0.5 is not negligible. It can be also seen that for 
£) = 4 we have /3i = 0, and (|189p starts with an term. 
The elapsed time as a function of the oscillaton mass can 
be obtained by expressing e from (|185l) and substituting 
into (ITM| . 
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For D = 3 spatial dimensions it is natural to start with 
a maximal mass configuration M — Mmaxj and wait for 
the mass to decrease until the ratio M /M^^.^ reaches a 
given value. Since the elapsed time t is inversely propor- 
tional to the scalar field mass m, in table IVlJ we list the 
product tm. 





£ 


tm 


t 

year 


eV 


0.01 


0.482 


5.35 


10i« 


1.12 ■ 


10"^ 


0.1 


0.383 


3.00 


10^' 


6.26 ■ 


10^2 


0.2 


0.320 


1.50 


1026 


3.12 


■ 10^ 


0.3 


0.269 


4.42 


10^' 


9.22 


■ 10** 


0.31884 


0.260 


6.57 


1032 


1.37 


10'« 


0.4 


0.224 


3.99 


1038 


8.32 


lO's 


0.5 


0.182 


1.22 


1048 


2.55 


1025 


0.6 


0.144 


1.28 


1062 


2.67 


1039 


0.7 


0.107 


1.94 


1085 


4.04 


1062 



TABLE VII: The time necessary for the oscillaton mass to 
decrease to M from the value Mmax a.tt = Q. The value of tm 
is given in Planck units, and also when the time is measured 
in years and the scalar mass in electron volts. 

Next we address the question that how much of its 
mass an initially maximal mass oscillaton loses during the 
age of the universe, which we take to be 1.37 • 10^" years. 
In Table I Villi we list the resulting oscillaton masses in 
units of solar masses (Mq), as a function of the scalar 
field mass in eV/c^ units. In order to facilitate com- 



eV 






M 


A/max -Af 
A/max 


lQ-35 


5.09 ■ 10" 


20 


8.20 • 10^'' 


1.41 ■ 10" 


38 


10-30 


5.09 ■ 10" 


15 


8.20 • 10'^ 


1.41 ■ 10" 


28 


10-25 


5.09 ■ 10" 


10 


8.20 • 10" 


1.41 ■ 10" 


18 


10-20 


5.08 ■ 10" 


-5 


8.20 • lO'-* 


1.40 • 10" 


-8 


10-'^ 


0.0704 




7.99 • 10* 


0.0258 




10-10 


0.163 




7.14 • 10"' 


0.129 




10"^ 


0.223 




6.30 • lO"'^ 


0.232 




1 


0.266 




5.58 • 10"'' 


0.319 




10^ 


0.297 




5.00 • 10""^ 


0.390 




10"' 


0.322 




4.52 • 10-2' 


0.449 




10"' 


0.342 




4.12 • 10^2'' 


0.498 





TABLE VIII: Mass M of an initially maximal mass oscillaton 
after a period corresponding to the age of the universe for 
various scalar field masses. The decrease in e from £max = 
0.525, and the relative mass change rate (A/max — M)/Mmax 
is also given. 

parison, we have chosen the same scalar field masses as 
in Eq. (178) of [3l|. The first two orders of the small- 
amplitude expansion yielded wil/ujax = 0.614 in Planck 



units for the maximal mass of the oscillaton. Taking 
the scalar mass in electron volts, this corresponds to 
Afmax = 8.20 • lO-"M0ey/(mc2), which was used in Ta- 
ble [VllTl The value mMmax — 0.607 from the numerical 
solution of the Fourier mode equations calculated in [51j 
corresponds to Afmax = 8.11 • lO^^^Mr^eV/ (toc^) in natu- 
ral units, which is the value used in l3l [. Comparing our 
Table rvml to the numbers in (178) of [3i|, after compen- 
sating for the shift in the initial mass, it is apparent, that 
for small scalar field masses, i.e. for m < 10~^^eV/c^ , os- 
cillatons decay more slowly in [3l|. The reason for this 
is that ^l[ uses a linear dependence of the mass on the 
small parameter, and consequently underestimates the 
radiation rate close to the maximum mass, similarly as 
we did in (I184p . For m > 10~^eV/c^ oscillatons radiate 
faster in ^31] , which is a consequence of the much larger 
value of the constant C3 in the mass loss law (|171[) used 
there. In spite of the differences, the overall picture re- 
mains essentially the same. For all scalar field masses 
that appear physically reasonable, a maximal mass oscil- 
laton loses a significant part of its mass during the life- 
time of the universe. This mass decrease is greater than 
10% if m > 4.57 • IQ-^^eV/c^, but it remains below 50% 
if m < 1.85 • lO^^eV/c^. The above results support the 
possibility that provided a scalar field exist in Nature, 
at least some of the dark matter content of our Universe 
would be in the form of oscillatons. 



VI. CONCLUSIONS 

We have derived an infinite set of radial ODEs de- 
termining the spatial field profiles of bounded solutions 
of time-dependent, spherically symmetric Einstein-scalar 
field equations in the limit when the scalar field am- 
plitude tends to zero. The lowest order equations are 
nothing but the D-dimensional generalization of the 
Schrodinger-Newton (SN) eqs. The SN eqs. admit glob- 
ally regular, exponentially decreasing solutions for spatial 
dimensions 2 < D < 6. The eqs. corresponding to higher 
orders in the expansion are linear inhomogenous ODEs. 
The class of solutions wc are interested in are oscilla- 
tons, which loose slowly their mass by scalar radiation. 
In the small-amplitude expansion we have obtained an 
asymptotic series for the spatially well localized core of 
oscillatons and related their radiation amplitude to that 
of the standing wave tail of exactly time-periodic quasi- 
breathers. For the class of symmetric scalar potentials 
we have determined the amplitude of the standing wave 
tail of time-periodic quasibreathers analytically adapting 
the method of Segur-Kruskal and using Borel summation. 
We have explicitly computed the mass loss rate for the 
Einstein-Klein-Gordon system in Z? = 3, 4, 5. 
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Appendix A: Einstein tensor and the wave equation 

The components of the Einstein tensor in the general spherically symmetric coordinate system (|12p are 



Gtt = 
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The wave equation (|4]) takes the form 



~B~ 



A 2AC"-^ 
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(A2) 
(A3) 



(A4) 
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(A6) 



Appendix B: Small-amplitude expansion in 
Schwarzschild coordinates 

In the main part of the paper we have used the spatially 
conformally flat coordinate system C = r^B. In this 
appendix we present the results of the e expansion in C = 
Schwarzschild area coordinates, in order to compare 
and to point out the disadvantages. The time dependence 
of the scalar field </> and the metric components A and B 
up to order are 



e^P2 COST + 0(6*), 



A = l+ e^a2 + e^af^ cos(2t) + 0{e^) , 



(Bl) 
(B2) 
(B3) 



where P2 , 02 , 02 ^^'^ ^2 are functions of p. The functions 
02 and P2 are again determined by the coupled differential 
equations and P5|) . resulting in the Schrodinger- 
Newton equations. However, &2 is determined as 



P 



D 



da2 
"2 'd^ 



(B4) 



instead of (|4T|) . The most important difference is the ap- 
pearance of the cos(2r) term in (jB2p . causing an or- 



der oscillation in the metric component gu- In spatially 
conformally flat coordinates there are only order os- 
cillating terms in the metric components. The amplitude 
of the oscillation is determined by the field equations as 



(2) 

a2 = -02 



(B5) 



Substituting into the expression (j72p of the magnitude 
of the acceleration of constant (r, 6*1, 02---) observers, to 
leading order we get 




cos(2r) 



(B6) 
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